Stretching Instability of Helical Springs 



David A. Kessler and Yitzhak Rabin 
Dept. of Physics, Bar-Ban University, Ramat-Gan, Israel 
(Dated: February 1, 2008) 

We show that when a gradually increasing tensile force is applied to the ends of a helical spring 
with sufficiently large ratios of radius to pitch and twist to bending rigidity, the end-to-end distance 
undergoes a sequence of discontinuous stretching transitions. Subsequent decrease of the force leads 
to step-like contraction and hysteresis is observed. For finite helices, the number of these transitions 
increases with the number of helical turns but only one stretching and one contraction instability 
survive in the limit of an infinite helix. We calculate the critical line that separates the region of 
parameters in which the deformation is continuous from that in which stretching instabilities occur, 
and propose experimental tests of our predictions. 

PACS numbers: 46.32.+X, 46.25.Cc,62.20.Fe 

The study of mechanical instability of thin elastic rods goes back to Euler and his contemporaries who considered 
the buckling of such rods under compression and torque j^, |2j. The present paper presents a theoretical prediction of 
a hitherto unknown instability that arises when a helical spring is stretched by a tensile force applied to its ends. As 
we will show in the following, stretching instabilities occur only in helices whose radius is sufficiently larger than their 
pitch and whose rigidity with respect to twist exceeds that with respect to bending. 

Consider a helical spring of contour length L, characterized by constant spontaneous curvature kq and torsion tq or, 
alternatively, by its radius r — ko/(ko + '^o ) pitch p — 2'ktq/{kq + Tq). Under the action of a constant stretching 
force F applied to its ends and directed along the z-axis (the ends are otherwise unconstrained and no torque is 
applied), the helix deforms into a curve whose shape is determined by minimizing the elastic energy]^, ^, |^, ^ 

^=^11°'/ ds[5u,{s)f ~ FR,, (1) 

with respect to the deviations 8L0i{s) of the generalized curvatures and torsions ijJi{s) from their spontaneous values 
in the undeformed state (s is the contour parameter). In general, these parameters are related to the curvature k, 
torsion t, and the angle a between one of the principal axes of the cross section (ti) and the binormal: cji = Kcosa, 
UJ2 = Ksina and 0^3 — t + da/ds. We study here the case that the only contribution to spontaneous twist comes from 
torsion (ao = 0) and therefore wqi = '^O: <^02 = and wqs = tq. The coefficients ai and 02 are bending rigidities 
associated with the principal axes of inertia of the (in general, noncircular) cross section and 03 is the twist rigidity. 
In the following we treat as given material parameters of the spring. The shape of the deformed spring can be 
obtained by inserting the {uji{s)} that minimize the elastic energy into the generalized Frenet equations for the unit 
vectors (t^ is the tangent to the curve at point s and ti and t2 point along the principal axes of symmetry of the 
cross section): 

dti/ds = - ^ gyfeWj-tfe, (2) 
jk 

where e^fc is the antisymmetric unit tensor0|. The space curve r(s) associated with the centerline of the deformed 
spring is then obtained by integrating the relation dv/ds = ts. The system of equations is closed by substituting the 
expression for the projection of the end-to-end vector on the z-axis, = dst^ls) ■ z, into Eq. |l|. Since the spring 
will always orient itself along the direction of the force, in the following we replace R^ by the end-to end-distance R. 

The force-extension curves are calculated numerically for a helix of 4 turns by finding a set ofuJi{sj) (the continuous 
curve is replaced by a discrete set of points, {sj} , j = 1,2,...) that minimizes the energy, Eq. |^. Two typical types of 
behavior are found. For sufficiently large ratio of pitch to radius p/r ^ 2'ktq/ and sufficiently small ratio of twist to 
bending rigidity as/ab (where a^^ = (aj~^ -I- a2^)/2), the deformation is continuous and identical curves are obtained 
by stretching the undeformed helix and by starting from a fully stretched configuration and decreasing the tensile 
force to zero (Fig. 1). 

In the opposite limit (sufficiently large r/p and as/at ratios), increasing F results in a sequence of discontinuous 
upward jumps of the end-to-end distance (Fig. 2), whose number depends on the number of helical turns and on the 
increment of the force. Comparison of the conformations on both sides of a discontinuity shows that the jumps are 
associated with the elimination of helical turns (Fig. 2, insets). Discontinuous downward jumps of the end-to-end 
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FIG. 1: Continuous deformation of a helix of 4 turns with ko = 1, to = 0.2, ai = 1, a2 — 5, as — 3. We show the end-to-end 
distance as a function of the force for both increasing (circles) and decreasing (solid hne) force. 



distance are also observed when one starts from a fully stretched configuration and gradually decreases the tensile 
force but the locations of these jumps are, in general, different from those of the upward ones. The observation 
of hysteresis loops of stretching and contraction suggests an analogy with first order phase transitions and implies 
that the elastic energy E{R) has multiple local minima whose depth and location vary with the tensile force F. In 
the absence of thermal fluctuations, a transition to a new energy minimum (not necessarily the lowest energy one) 
takes place when the minimum corresponding to the original state disappears. While in systems undergoing phase 
transitions the choice of the final state is determined by the kinetics of the temperature quench, in our case it depends 
on the magnitude of the increment of the force. 

The analogy with phase transitions suggests that there exists a critical surface in the parameter space that separates 
between regions in which the deformation is continuous from those in which a discontinuous stretching transition takes 
place. It is of particular interest to locate this critical surface in the limit of an infinitely long helix, for which end effects 
are negligible and it is expected that the only stable states are perturbations of the original undeformed helix and of 
the completely stretched spring (and therefore only one hysteresis loop that corresponds to a single stretching and a 
single contraction instabilities, is expected). Insight into this limit can be achieved via analytical approximations. In 
order to reduce the number of variables, we consider a ribbon with highly asymmetric cross section, a-ijax —>■ oo, in 
which case the cross section of the ribbon is pinned to the Frenct frame of its centerline, a{s) — 0, and the deformed 
spring is completely characterized by k and r. Based on our numerical studies of short helices, we would like to stress 
that varying 02 (while keeping all other parameters fixed) affects the detailed extension versus force curves but did 
not change the character of the transitions, and stretching instabilities accompanied by hysteresis were observed even 
in the symmetric case, ai — a2. 

In the following we measure the energy in units of uiKqL, and force in units of ciKq. In order to obtain analytical 
results, we consider the limit of an infinitely long helix and assume that the spring maintains a helical shape and 
responds to deformation only by adjusting its curvature and torsion parameters. Thus, we look for constant (s- 
independent) parameters k and r that minimize the energy per unit length, 

j^iT 1 / / 1^2 I 03/01 2 Ft/kq 

E/L = - [K/Ko - 1) + {t/kq - To/Ko) , ^ ^ - (3) 

Support for this variational approximation comes from numerical minimization of the energy, Eq. |lj which shows 
that while the curvature and torsion oscillate about their mean values, the amplitude of oscillation decreases with the 
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FIG. 2: Hysteresis loop for a helix of 4 turns with kq = 1, to = 0.2, a\ 
as a function of the force for both increasing and decreasing force. 



1, 02 — 5, a-i — 4. We show the end-to-end distance 



length of the spring. In Fig. 3 we compare the extension vs. force curve for a heUx of 12 turns, obtained by exact 
numerical minimization of Eq. ^ with that obtained using the variational estimate. The fact that the two curves 
coincide up to the stability limits of the initial and the final states provides further support for the validity of our 
variational approach. 

Analytical results can be obtained in several limiting cases. In the case of large twist rigidity, a^/ai — > oo, 
the torsion is pinned to its spontaneous value (r = tq) and the curvature obeys the relation ai {k/kq — 1) -I- 

r 2 2i"^/2 

F (to/kq) (i^/kq) (k/kq) + {t/kq) — 0. Graphical analysis of this equation shows that it admits either a single 

minimum or two minima separated by a maximum. For tq > Tq^c there is a single minimum of k as a function of F. 
For Tq < To,c, there is a window of F between which there are 3 solutions and outside which there is only one solution. 

5 /3 

This window closes at the critical point: tq^c/kq = (2/3) = 0.363, Fc = 0.651, and Kc/kq — 0.444. Another limit 
in which the critical point can be calculated analytically for arbitrary a^/ai is tq/kq 1. In this case we find that 
two minima exist for 03/01 > 4/3. The critical value of this ratio is therefore 4/3 and at this point F^ = 4/3, and 
Kc/kq — 54^/'°^ (to/kq)^^^ /2. In between these two limits the critical points have to be found numerically, resulting in 
a critical line in the 03/01 — Tq/kq plane (Fig. 4). 

We have shown that there exists a range of spontaneous curvatures and torsions and bending and twist rigidities 
in which a helical spring does not deform continuously with the tensile force. Instead, as the force is increased, the 
spring undergoes a sequence of stretching instabilities. Hysteresis is predicted to take place when the force is decreased 
starting from full extension, as both the number of sudden contractions and their locations differ from that of the 
stretching instabilities at increasing force. We presented a variational calculation of the critical line that separates 
the region of parameters in which the deformation varies continuously with the force from that in which stretching 
instabilities take place. Even though the analysis was restricted to the limit of an infinitely long spring with a highly 
asymmetric cross section, our numerical results indicate that stretching instabilities will be observed in helical springs 
with both the radius to pitch ratio and the ratio of twist rigidity to bending rigidity larger than some numbers of order 
unity. While the former requirement can be realized by preparing a spring of a desired shape, the latter condition 
may require the use of materials with special elastic properties. Notice that for isotropic materials the rigidities can 
be expressed in terms of the shear and Young moduli fi and E and the geometry of the cross section; for example, in 
the case of rectangular cross section with sides ci and C2 one has|^ oi = £'0102/12, 02 = EC1C2/I2, and 03 ~ /iCiC2/3 
(the last expression holds in the limit ci S> C2). For incompressible materials this yields 03/01 ~ 4/3 which coincides 



4 



80 



CD 



o 



40 



20 




— Constant K Model 

- - Increasing Force 

Decreasing Force 



0.5 



1.5 



Force 



FIG. 3: Hysteresis loop for a helix of 12 turns, with kq = 1, to — 0.2, ai = 1, a2 = as — oo. We show the end-to-end distance 
as a function of the force for both increasing (dashed line) and decreasing (dotted line) force. Also shown is the result of our 
simple model with constant k (solid line). 
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FIG. 4: Critical line in the as/ai — to/kq plane for an infinitely long helix with a2/ai — > oo. The vertical and horizontal dotted 
lines correspond to the two limits discussed in the text, ro,c//«o = 0.363 and 03/ai = 1.333, respectively. 
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with the critical value of the ratio. The ratio 03 / ai can be increased beyond this critical value by reducing the Poisson 
ratio or by using anisotropic materials with high resistance to twist. 

We would like to comment on the connection between the stretching instabilities of helical springs and those familiar 
from other branches of physics. The occurrence of multiple energy minima and hysteresis bring to mind first order 
phase transitions but since thermal fluctuations and entropy play no role in our work, the analogy is somewhat 
misleading. The predicted transitions appear to be more closely related to hydrodynamic instabilities, even though 
the latter arise due to nonlinear effects, while the former can be described within the framework of the linear theory 
of elasticity!^. Further reflection shows that the elastic theory of thin rods is truly linear only in the "Lagrangian" 
frame of the deformed object (i.e., in its intrinsic coordinate description). In the presence of external force one must 
transform into the laboratory frame (e.g., Euler angles) and the energy is no longer a quadratic form in the deviations 
from the undeformed state (recall that hydrodynamics is also a linear theory in the Lagrangian frame!). 

Finally, we would like to stress that since the present analysis is purely mechanical, the results are directly applicable 
only to macroscopic objects. The relevance of the present work to the deformation of nano-objects such as biopolymers 
(DNA, proteins), protein filaments (actin, microtubules, etc.), and carbon nanotubes depends on the question of 
whether these objects possess spontaneous curvature and torsion. Interestingly, discontinuous stretching transitions 
were reported in single molecule extension studies of torsionally constrained DNA at high degree of supercoiling 
(unlike our model, supercoiling was not spontaneous but was achieved by the application of torque) ||l^ . The study 
of stretching instabilities of such objects requires consideration of thermal fiuctuations and extension of the present 
analysis using the methods of references [Q and which we hope to report on soon. 
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